Optimization of the residence time of a Brownian particle in a 

spherical subdomain 

0. Benichou and R. Voituriez 

Laboratoire de Physique Theorique de la Matiere Condensee (UMR 7600), 
O ' case courrier 121, Universite Paris 6, 

4 Place Jussieu, 75255 Paris Cedex 

I— s 
^1" 



ctf 



O 

o 



X 



(Dated: July 6, 2010) 

Abstract 
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O ' In this communication, we show that the residence time of a Brownian particle, defined as the 

| j c— ,e tim e spe lrt ta a gi ve, regi oa o f space, can b e op tim « a S a tocta o f th e diffusi oa 

coefficient. We discuss the relevance of this effect to several schematic experimental situations, 
■*-» . 

classified in the nature - random or deterministic - both of the observation time and of the starting 
position of the Brownian particle. 
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Among the numerous features of stochastic processes, the time it takes a random walker 
to reach a target - the so-called first-passage time (FPT) - is a crucial quantity that governs a 
variety of physical systems^"-. Indeed, numerous real situations, ranging from (sub) diffusion 
limited reactions^ to animals searching for foocr^ can be rephrased as first-passage prob- 
lems. In all these situations, the FPT is a limiting quantity, whose optimization is a crucial 
issue. 

Very often however, the reaction process is not infinitely efficient, so that the relevant 
quantity which has to be optimized is not the FPT to the target itself, but rather the 
residence time of the random walker in the vicinity of the target 14 . The importance of this 
observable, defined as the cumulative time spent by the random walker in a given subdomain 
centered on the target up to a fixed observation time (see figure [T]), comes from the fact 
that it can be seen as a measure of the interaction time between the random walker and the 
target. The study of the statistics of this general quantity has been a subject of interest for 
long, both for mathematicians 1 ^— and physicists^— . As a matter of fact, the residence 
time has proven to be a key quantity in various fields, ranging from astrophysics 25 , transport 
in porous medial and diffusion limited reactions 1 ^ 1 ^ 2 ^. 

The question we address here is to determine how the residence time depends on transport 
properties of the random walker, and more precisely on its diffusion coefficient. As we 
proceed to show, this dependence is non trivial, and in some situations it is possible to tune 
the value of the diffusion coefficient in order to maximize the residence time. 

To get an intuition of this effect, let us start with a simple analysis of the problem. 
Each time the random walker enters the interaction zone of typical size r, it spends inside 
a typical time of order r 2 /D which is a decreasing function of the diffusion coefficient D. 
On the other hand, after an observation time t, the number of times the interaction zone 
has been visited is an increasing function of D, and behaves like a/ 'Dt/r 2 in dimension 1, 
ln(Dt/r 2 ) in dimension 2 and like a constant in dimension d > 3, in the long time limit 
t — y oo^ 2 -. For a walker starting from the exterior of the interaction zone, it is thus clear 
that both limits of small and large diffusion coefficient D lead to a small residence time, 
indicating the existence of a maximum as a function of D. However, it has to be noted that 
the observation time has to be finite to make this optimization possible. Indeed, the mean 
residence time [i goes to infinity with the observation time in dimension d < 2, and tends 
to a constant divided by D in dimension d > 3, which is a decreasing function of D. 



With the exception of 28 - which only deals with a one-dimensional situation in the context 
of DNA/protein interactions - and o&tM - where it has been proposed that slow diffusion 
could be a way of enhancing the kinetics of imperfect reactions - the dependance of \x with 
D does not seem to have received a lot of attention so far, probably for two reasons. First, 
in most of calculations of residence times, the diffusion coefficient is taken fixed (D = 1/2 in 
the mathematical literature, D = 1 in^). Second, the observation time is generally taken to 
be infinite or at least large in explicit determinations of residence times (see for instance the 
Section III of 34 ), which, as mentioned above, does not permit to reveal the non monotonic 
behavior we are looking for. 

In the following, we focus on the generic case of a Brownian particle evolving in a 3- 
dimensional space, and discuss three "experimental" situations, potentially relevant to dif- 
ferent problems of chemical reactivity : (i) the observation time is a deterministic variable 
; (ii) the observation time is a random variable ; (iii) both the observation time and the 
starting position of the particle are random. 

Denoting R the initial position of the particle and r the radius of the target, assumed to 
be a sphere centered at the origin O, we study the pdf Pt(T) of the residence time T after 
an observation time t. Following Kao 1 ^, we introduce the following trapping problem : we 
assume that the Brownian particle disappears on the target with a constant and uniform 
rate k. The survival probability S(t) after an observation time t can be written in two 
different ways. First, it reads 

S(t)= f P t (t')e- kt, dt', (1) 
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where e~ kt represents the survival probability at time t' conditioned by the fact that the 
particle has spent a time t' is the reactive zone. Second, it is easy to see that S satisfies the 
backward Fokker-Planck equation 2 : 

d t S(t) = [DAr - kH(r - R)]S(t), (2) 

where H(x) = 1 if x > and otherwise, the laplacian Ar, involves derivatives with respect 
to the starting point R and R = |R|. This partial differential equation has to be completed 
by the initial condition S(t) — > 1 if t — > and the boundary condition S(t) —$■ 1 if R — ¥ oo. 
Denoting by S(s) = f^° e~ st S(t)dt the Laplace transform of the survival probability and 
Laplace transforming Eq. (J5J), it is easily seen that, if the starting point R is exterior to the 



sphere of radius r, 



while if R is interior 



S(s) = - s + AR x - d l 2 K dl2 _ x (rJ^\ , 
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where the constants A and B are given by continuity conditions of the survival probability 
and of its first spatial derivative at R = r, K v and I v stand for modified Bessel functions and 
the space dimension is d. Note that, in Eqs.03]),(JI|), only one of the two Bessel functions Id/2-1 
and Kd/2-1 generating the set of solutions appears, in order to fulfill boundary conditions 
when R — > and R — > 00. In our case d = 3, it is explicitly found that, if the starting point 
is exterior to the sphere of radius r, 



S(s) 
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k sinh (r y/(s + k)/D) - r y/(s + k)/D cosh (r y/(s + k)/D) 

s{s + k) ^JJjD sinh(r v /(s + k)/D) + ^(s + k)/D cosh(r yj(s + fc)/£>) 



(5) 



while if it is interior: 

1 sinh(iV(s + k )/ D ) 
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k(l + r^/sjD) 
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p(p + k ) y/s/Dswh(ry/(s + k)/D) + y/[s + fc)/£> cosh(r^(s + k)/D) 



(6) 



The Laplace transform j2(s) of the mean residence time /x(t) is then obtained by expanding 
Eqs(jni)-(ISD as a function of the parameter k, since S(s) writes: 



S(s) 



dte~ st / &t'e~ kt P t (t') = -- kj2(s) + o(k) 



(7) 



If the starting point is exterior to the sphere of radius r, the small k expansion of EqfJSJ) 
leads to : 



fi(s) 
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' 's/D cosh(r \/s/D) — sinh (ry/s/D) 



(8) 



while if the starting point is interior to the sphere of radius r, the small k expansion of EqflBJ 
leads to : 



1 (^ (l + r^/s~/D)smh(R^/s~/D)' 
Ry/s/Dexp(ry/sJB) 



M s ) = 72 I l » /-m — ; — ^^ — I ' ( 9 ) 



Deterministic observation time. In the first case of a deterministic observation time t, 
corresponding for instance to a chemical reaction designed to be terminated after a time t, 
it is actually possible to Laplace invert the transform of the mean residence time. Using the 
two known Laplace inverts: 

C-\e-^/s-){t) = -^-e-^ + l - (l-erf(^)(2t+l) (10) 

£ -i (e -Vi /s 5/ 2)(0 = _1 ^ _ erf (J_yj (1 + 6t) + 1 yie-V(^) (4 t + 1) (11) 

where £ _1 stands for the inverse Laplace operator and erf (x) is the error function, a direct 
Laplace inversion of Eqs.flU),© gives the explicit expression 
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Note that in the limit of large observation times t — > oo, this expression gives back the known 

results of— : w(oo) = — — if R < r (compare with Eq.(3.16a) of 34 ) and /i(oo) = 

6D 3RD 

if R > r (compare with Eq.(3.11) of^ 4 -). Moreover, as expected, it is easily seen that, if 
R > r, /i(t) — y if D — y or if D — y oo. As a consequence, there exists a value of Dq that 
optimizes the residence time (see figure [2]). Deriving the expressions of n(t) with respect to 
D, this optimal value is easily seen to satisfy the following implicit equation: 

_ 2 e " *to$~ v /l^i (# 2 _ 2 r 2 - i?r - 2 D t) + 2 e -1 ^^ y^t (# 2 + R r - 2 D t - 2 r 2 ) 
-v^erf^^^^(i?-2r)(i? + r) 2 + v ^erf^^±|^)(2r + J R)(-r + J R) 2 = (13) 



whose solution can be discussed in limiting regimes. When R ^> r, we find that D Q ~ aR 2 /t, 
where a ~ 0.353 is the implicit solution of the equation 



erf \-X= 1 - 1 



p -l/(4a) 

+ — -0, (14) 
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obtained by replacing D by aR 2 /t in Eg. (1131) and taking the limit R^> r. In other words, 
the optimal diffusion coefficient D does not depend of the extension of the target r in 
this regime. This is in strong contrast with the opposite regime R — >■ r + , which leads 
after a similar analysis to an optimal diffusion coefficient D$ ~ 3r 2 (R/r — l)/(2£), which, 
in particular, vanishes when i? — >■ r + as could be expected from the qualitative analysis 
presented above. 

Random observation time. We now consider another experimental situation, correspond- 
ing to a random observation time, distributed according to an exponential law with mean 
1/p. This case corresponds for example to a situation where the particle has a finite lifetime, 
as in many biological situations. Averaging over this random lifetime, the mean residence 
time becomes: 



(A*) 



/*oo rt 

/ dtpe~ pt / dt't'P t (t') = P j2(s = p), (15) 
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where ju(s) is given by Eqs.(J21)-(!EJ)- Once again, for R > r, the mean residence time can 
be shown to have an optimum as a function of D (see figure [2]). Interestingly, the optimal 
diffusion coefficient D satisfies the simple implicit equation 

dx -\- x 

tanhx = — , (16) 

x l + ax + 1 

where x = r\fpj ' Dq and a = R/r. 

For a ^> 1 , the Lagrange inversion formula^ allows one to obtain explicitly the expansion 

of x as powers of 1/a, up to arbitrary order. The first terms are given by: x = ^ + -^ + .... 

In the opposite limit a — > 1, it is easy to deduce from Eq. (fT6|) that x ~ . Finally, 

\Ja — 1 
this gives again two regimes for the optimal diffusion coefficient: Dq ~ pR 2 /A when R^> r 

and D ~ pr 2 (R/r — 1) when R — y r + , in close analogy with the two regimes found in the 

case of a deterministic observation time. 

Random observation time and random initial position. In the last situation, we assume 

that both the observation time and the initial position are random variables. More precisely, 

the observation time is again assumed to be distributed according to an exponential law with 



mean 1/p, while the initial position is assumed to be uniformly distributed in an annulus 
region centered on 0, of inner radius b and outer radius c. This mimics for example a case 
of chemical reactions where the position of the source of reactants cannot be accurately 
determined. In this case of an annulus region, it is not clear a priori if the mean residence 
time is still non trivially optimizable with respect to D. Indeed, if the initial position is 
within the reactive sphere, the best diffusion coefficient is simply Dq = 0. On the contrary, 
if the annulus region is entirely exterior to the reactive sphere, a non trivial optimal Dq is 
expected from previous results. This raises the question of the intermediate case where the 
annulus partially covers the reactive sphere. 

In the first case of an annulus region exterior to the target zone (b > r), we find 



W) M 



sinh(r/Z) ' 

- cosh(r/i) 



p(c 3 — 6 3 ) r/l 

x [( c + l) e - c / l -(b + l)e- b/l ], (17) 

where / = \J D/p. This quantity is easily shown to tend towards when D — > oo while 

This proves that (/x) admits an optimum diffusion coefficient, as expected in view of the 
previous results. 

In the opposite situation of a annulus region such that < b < r < c, we get 

1 r 3 - 6 3 _ 3/ 

W ~ p-^Ttf ~ p( c 3 _ 53) X 

[e-' r/l (r + /) (b cosh (b/l) - I sinh (b/l)) 

-e- c/l (c + /) (r cosh (r/l) - I sinh (r/l))] . (19) 

The local behavior of this quantity when D — > actually depends on the relative order of 
the two distances r — b and c — r. If r — b < c — r, 

7JA _ l - r lzl_ Zrbl {r _ h)/l 

W pc 3 -6 3 2p(c 3 -6 3 ) ' l } 

showing that there is an optimal diffusion coefficient. 
Remarkably, ifr — b > c — r, 

M - -^ ~ -T^\/^T e " (c " r)v/ ^ < °> ( 21 ) 

pc 6 — b 6 p(c d — o J V p 2 




and this time there is no possible optimization. 

The last case of an annulus region of inner radius equal to zero (i.e. b = 0), i.e. a 
sphere of radius c, deserves some specific attention. In this case, Eq. (FIT)]) is still valid but, 
quite surprisingly, no optimization with respect to the diffusion coefficient is possible. This 
monotonic decrease is compatible with the local behavior (D — > 0): 

V 

We stress that this result holds for arbitrary values of the radius c. 

In conclusion, we have proposed a very simple mechanism that allows one to optimize 
the residence time of diffusing molecules with respect to the diffusion coefficient. We have 
discussed the relevance of this effect to several schematic experimental situations, classified 
in the nature - random or deterministic - both of the observation time and of the starting 
position. We believe that this optimization of the residence time discussed in the case of 
the standard Brownian motion is robust (see in particular Supplementary Material for a 
discussion of a discrete space version of the model presented here) and is generalizable to 
more complex transport processes. 
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Figure captions: 

Figl: A Brownian particle crossing an interaction sphere. The residence time is the 
cumulative time spent by the particle inside the sphere. 

Fig2: The mean residence time as a function of the diffusion coefficient for deterministic 
(t = 1) and random observation times (p = 1). In both cases, the Brownian particle starts 
outside of the sphere (R = 3, r = 1). The optimal diffusion coefficient in the deterministic 
case (resp. random) is well approximated by the limiting expression Dq ps aR 2 /t (resp. 
D pa pR 2 /4) given in the text, even if here the condition of applicability R^> r is not well 
satisfied. 
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FIG. 1: A Brownian particle crossing an interaction sphere. The residence time is the cumulative 
time spent by the particle inside the sphere. 
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FIG. 2: The mean residence time as a function of the diffusion coefficient for deterministic (t = 1) 
and random observation times (p = 1). In both cases, the Brownian particle starts outside of the 
sphere (R = 3, r = 1). The optimal diffusion coefficient in the deterministic case (resp. random) 
is well approximated by the limiting expression Dq ~ aR 2 /t (resp. Dq ~ pR 2 /4) given in the text, 
even if here the condition of applicability R 3> r is not well satisfied. 
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